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We propose a set of interferometric methods on how to detect Majorana bound states induced by 
a topological insulator. The existence of these states can be easily determined by the conductance 
oscillations as function of magnetic flux and/or electric voltage. We study the system in the presence 
and absence of Majorana bound states and observe strikingly different behaviors. Importantly, we 
show that the presence of coupled Majorana bound states can induce a persistent current in absence 
of any external magnetic field. 
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I. INTRODUCTION 



One of the main interests in current research on quantum computation is to find new materials that facilitate the 
'~~| ■ construction of large scale quantum computers. The main impediment is decoherence and the generation of errors. A 
d \ promising way of reducing the effect of decoherence and errors is to employ quantum systems that have topological 
' characteristics^. Recently, topological insulators have been considered that can be tuned to support topological 
lyj \ states, such as Majorana bound states. There have been several proposala^i^ how one could encode information with 
D . Majorana fermions that is protected against a variety of errors. Majorana fermions which can also occur in highly 
^ ■ correlated systems like Px + ipy wave superconductors^, the v = b/2 fractional quantum hall stated, at the boundary 
I ; , of supcrfluid 3Hc-BX and finally in superconducting graphene^ have the special characteristic that they are their own 
■ antiparticles. They are also predicted to appear as low energy excitations in Kitacv's two-dimensional spin-1/2 system 
^ ' on a honeycomb lattice^i^i^. 

I ' The aim of this work is to detect Majorana bound states, at the interface between topological insulators with 
\ superconducting and magnetic correlations, addressed also in Refs.(6lll0j. Coupled Majorana bound states (implying, 
. two Majorana bound states which are interacting) can encode a qubit non-locally and obviate local environmental 
Q ' pcrturbationsiS. Their detection, so far, has been difficult as Majorana fermions are neutral quasi-particlcs. To 
overcome that recent works proposed to employ Dirac to Majorana fcrmion converter o^^'^^ . There, Dirac particles 
emitted from a source are converted to Majorana's and then reconverted back at the drain. In our work Majorana 
bound states are efficiently monitored by a mesoscopic Aharonov-Bohm interferometer. In particular, the presence 
of Majorana bound states can be probed by the symmetry of the non-local conductance as function of the magnetic 
' field or an applied electric field. Further, we show that the presence of coupled Majorana bound states could induce 
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FIG. 1: (Color online) An overview of the setting from the top. The 2D topological insulator (TI, in blue) is made into 
an Aharonov-Bohm interferometer (a). In the upper arm of the ring, MBS (white ellipses at edge of Ferromagnet(F) and 
Superconductor(S)) occur as shown in (c). Topological edge modes (EM), circulate along the edges, interacting with the MBS. 
A magnetic flux penetrates the ring while an electric voltage covers the lower arm (in green) of the ring, (b) A loop made up 
of a TI which supports MBS. (d) Representation of incoming and outgoing waves from left coupler as in (a). 
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persistent currents in a topological insulator ring in absence of any magnetic flux. 



II. THE PROPOSED PHYSICAL MODEL 



To detect Majorana bound states(MBS) that exist in the topological insulator(TI)^, we look at an Aharonov-Bohm 
(AB) interferometer made up of a 2D TI (e.g., Bi2Se-j^ or a HgTc quantum wclli^) as depicted in Fig. 1(a). In a TI, 
spinorbit coupling causes an insulating material to acquire protected edge or surface states. A magnetic flux exists in 
the center of the AB interferometer. The regions of the interferometer arc labeled in Fig. 1(a). In the upper arm of 
the ring, at the interface between s-wave superconducting(S) and a thin ferromagnetic layer(F), MBS (white ellipses 
at the interface) occur as shown in (c). Topological edge modes(EM), circulate at the edges, interacting with the 
MBS. A magnetic flux penetrates the ring while an electric voltage covers the lower arm (in green) of the ring. The 
superconductor and Ferromagnet deposited on top of the TI, via the proximity effect, induce superconducting and 
ferromagnetic correlations in the substrate. The places where these correlations intersect is where MBS occur. The 
ring is connected to two leads on either side. The left lead is at potential Vi while V2 = 0(see Fig. 2). The full 
Hamiltonian in the upper and lower arms of the interferometer satisfies 

{vpr^cr^ + [eV - Ep + eA/hc)T^)^ = , (1) 

wherein p = —iTid/dx is the momentum operator, Ep the Fermi energy, v the Fermi velocity, eV the electric field 
applied to the lower arm only and A defines the magnetic vector potential. The four component wave-function 
^' — (^fg^, ^'ej., 5'/4)"^, while the t matrices mix the e and h blocks of the Hamiltonian. The eigenstates of 
Hamiltonian ([T]), can be calculated by considering plane wave solutions. The Hamiltonian for the superconducting- 
magnet interface (white ellipse at edge of red and greens dashed areas as in Fig. 1(c)) is that for the MBS, 

Hm = —(jyEm (2) 

. As discovered by Fu and Kane»^, a MB state appears at the intersection of the magnet-superconductor interface 
with the edge of the TI. The 4x4 S-matrix of scattering via the MBS, Ref. [l^l, can be written as SMaj = sfj where 
{a, b} = {e, h} and {i,j} = {1, 2}. The elements are determined as follows 
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In the above equation, E is the incident electron energy, ri/2 are the strengths of coupling to left/right arms and Em 
the strength of coupling between the individual MBS. 



III. EDGE MODES AND SCATTERING MATRICES 



To understand the edge modes and the type of scattering matrix needed to describe them, we show in Fig. 2, the 
edge modes flowing in our system and the way MBS affect them. We start with a related system, in Fig. 2(A) a 
quantum hall conductor with an AB flux is shown. A localized state around the hole develops and is sensitive to the 
flux, the outer edge states are nol^i^. However the outer edge states determine the net conductance when a voltage is 
applied to left lead. Thus if there is no scattering between outer and inner edge states the conductance is independent 
of flux. For conductance to be sensitive to flux one has to couple both these states. Thus as seen in Fig. 2(B) we 
couple the outer edge states to the inner states via the two couplers. The couplers induce inter edge scattering as 
shown by dashed lines. A 3 x 3 matrix can effectively describe this scattering process as there are 3 outgoing and 3 
incoming modes. In the TI with a hole shown in Fig. 2(C), the edge states occur in pairs but are of opposite spins 
(black:up, white:down), there are also two localized counter-propagating edge states which develop around the hole. If 
as we assume there is no scattering between states with opposite spins then up spin conductance will be equal to that 
for down spins. Couplers induce backscattcring for all edge modes while MBS mixes electron and hole edge modes and 
also backscattcrs. We restrict ourselves only to describing electron spin up and hole spin down edge modes as shown 
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FIG. 2: (Color online) Edge modes and back scattering. (A) Quantum Hall conductor with localized flux. (B) Coupling 
between outer and inner edge modes in a Quantum Hall conductor. (C) A TI with a hole, spin -up(black) and down(white) 
edge modes counter propagate along the hole. (D) A MBS scatterer (red shaded area on the upper arm of the AB ring) induces 
e-h mixing and back scattering. 




FIG. 3: (Color online) Gni, in units of e'^ /h, as function of incident electron energy E (controlled via a gate voltage) in units 
of A. The dashed and solid lines in the left panel are for individual {Em = 0) and coupled MBS {Em = 1-0), Em in units of 
A the superconducting gap. 

in Fig. 2(D) this leads to two counter propagating electron and another two hole edge modes around the flux. This 
leads to a 6 X 6 S-matrix for either couplers while a 4 x 4 matrix couples inner edge modes. Similarly one can construct 
the electron spin down and hole spin up edge modes. The total conductance will be twice that for incident electron 
spin-up edge modes since there is no spin-flip scattering. In Fig. 1(a) and 2(D), the length of upper arm is lu and of 
lower arm is Id, total circumference of loop is L = -|- Id- The Majorana scatterer further divides the upper arm, as 
shown in the figure /„ ^h+h with li ~ h- The loop is connected to two current leads on either side. The couplers 
(triangles) in Fig. 1 which connect the leads and the loop are described by a scattering matrix S. The S matrix for the 
left coupler yields the amplitudes Oi = {Se, Sh,"fie,"/ih, Pdei Pdh) emanating from the coupler in terms of the incident 
waves h = {oe,Oh,ade,adh,Pie,Pih), and for the right coupler yields the amplitudes O2 = {Te-,Th.,^de-,ldh, I32e, P2h) 
emanating from the coupler in terms of the incident waves I2 = (*ej ^de, Cd/u 72e, 72h)- The S-matrix for either of 
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the couplersi^, left and right, is given by 

/-(a + 6)1 yel yel\ 
S = \ al bl \ (4) 

\ yil bl al J 

with a = (1 — 2e) — = — 2e) + 1), I being the identity matrix, e plays the role of a coupler with 

maximum coupling ^ — \ while for e = 0, the coupler completely disconnects the loop from the lead. For left coupler, 
the waves into and out are marked in Fig. 1(d). The waves incident into the branches of the loop are related by the 
S matrices for left part of the upper branch by 
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with k'^ = {E + Ef + V)/'hv, k'f^ = {E — Ef — V)/fiv are the electron and hole wave- vectors, while fcg and kh are 
wave- vectors with = 0. ^ and ^ are the phase shifts due to flux in the upper and lower branches. Clearly, 
^ -t- ^ + ^ = where $ is the flux piercing the loop and $o is the flux quantum The transmission 

and reflection probabilities from Fig. 1 and Eq. (3) are given as follows: normal electron reflection R^. = non- 
local electron co-tunneling Tg = local Andreev reflection = and non-local crossed Andreev reflection 
Th = wherein t,6 are as depicted in Fig. 1(a). In the calculations we consider e = h = c = 1 and Fi = r2 = F. 
For the setting as described in Fig. 1(b), the scattering matrices can be written in exactly similar fashion (see also 
Ref.[l^). The total persistent current densitj*ii in a small interval dE is then sum of the individual electronic and 
hole current densities calculated as J = ev{Je + Jh), with Je = (|7ieP — l/^ieP), and Jh = (I71/1P — |/3i/ip). 



IV. DETECTION PROCEDURE FOR MBS 



To distinguish the behavior of the MBS we analytically solve Eqs. (4-7) and derive expressions for reflection and 
transmission probabilities. For brevity we present plots that result from these solutions. In Figs. 3 and 4 we plot the 
nonlocal conductance G„; = {e'^ /2h)[l- Re- Rh+Te + Th] = {e^ /h){Te+Th). We subsumed a "-" sign into the hole 
wavevectors, indicating their opposite direction to electrons, in the S matrices. Thus in Gni we add the individual 
contributions rather than subtract. As a consistency check- the sum of probabilities T^ + T^ + Re + Rh — 1- Further, 
current conservation also holds as currents in either leads are equal T^+Th = I — Re — Rh- The non-local conductance 
implies a current which appears in the right lead, while a voltage Vi is applied to left lead and no voltage is applied to 
the right. Hence, appearance of a current is due to the non-local effect of voltage applied to the left lead. In Figs. 3-5, 
the dimcnsionless parameter e 1/2, Fermi energy Ep — and F = 0.1 in units of A and lengths 1^ = 1^^ = 1/2 in 
units of hv/A. In Fig. 3, we focus on the behavior of G„; as function of the electronic energy (which can be tuned 
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FIG. 4: (Color online) Non-local conductance, in units of /Ti, as function of magnetic flux (top panel) in units of $o = fi-c/e 
and electric potential (bottom panel) in units of A for E = 0.1 A. 
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FIG. 5: (Color online) The current density(J) in units of ev as function of the incident electron energy, for e = 1/2, and Ef ~ 0. 
A finite current (red dashed line) flows when MBS are coupled. 



by a gate voltage) in the top left panel. We see there is a pronounced dip in the case wherein MBS are decoupled. 
On the adjacent panel we plot the individual contributions. We see that most of the contribution to Gni comes from 
electron co-tunneling. Further, Gni{E) = Gni{—E) for = 0, but this equality does not hold in presence of magnetic 
flux. Fig. 4 shows the variation of G„j as function of the magnetic flux. It shows Gni{4>) = Gni{—(t>) in the case 
where the MBS are decoupled, while for coupled states Gniifj)) 0). The absence or presence of an electric 

voltage on the lower arm does not make much of a qualitative difference while a quantitative difference is manifest. 
On the lower panel of Fig. 4 we plot Gni against the electric potential V on the lower arm of the ring. It is seen that 
while Gni{V) ^ Gni{—V) irrespective of whether MBS are coupled or not, there is a halving of periodicity when MBS 
are decoupled. The reason for non-observance of magnetic field symmetry is two-fold: (i) A coupled MBS scatterer. 
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FIG. 6: (Color online) Non-local conductance as function of magnetic flux and voltage applied to left lead, Vi, for a AB ring 
with a s-wave superconductor in its upper arm in absence of MBS. z represents normal metal-superconductor interface strengths 
in units of A. 



TABLE I: Detecting Majorana bound states (MBS) 



MBS 


Magnetic field 


Electric field/ gate voltage 


individual {Em = 0) 


G{<f>) = G(-0) 


G{E) = G{^E) 


coupled {Em / 0) 


G(<^) / G(-0) 


G{E) = G{-E) 


Absent 


G(</>) = G{-4>) 


G{Vi)^G{-Vi) 



breaks time reversal symmetry as in Eq. e.g., sf| ^ and (ii) breakdown of Andreev reflection symmetry 
{sli*{—E) = Sif{E)) instead of sli*{—E) = —Sif{E), regardless of whether MBS are coupled or not. However, 
for E,r A, the breakdown of Andreev reflection symmetry can be discounted as this amplitude is minimal and 
therefore we chose the low energy and weak coupling sectors in Figs. 3-4 so as to bring out the distinction between 
coupled and individual MBS in the non-local conductance versus magnetic field plots wherein breakdown of time 
reversal symmetry is the main reason. 

Because, of break down of time reversal symmetry electrons and holes scattered from the coupled state get different 
phases when they travel from left to right and vice-versa. The same philosophy which first predicted persistent 
currents in mesoscopic rings due to the fact that a magnetic flux can break time reversal symmetr}*i^ is once again 
present here. Importantly, here it is not just the magnetic flux which causes the breaking of time reversal symmetry 
but it is also the Majorana scattering. This implies that a circulating current can arise in the TI loop because of 
the scattering due to coupled MBS independent of the fact whether a magnetic flux is present or absent. To isolate 
this effect we calculate the persistent current for the setting described in Fig. 1(b). In presence of coupled MBS, a 
persistent current is induced in the TI loop while for individual (or. de-coupled) MBS such a current is absent. In 
Fig. 5, we plot the persistent current density (in units of eu), which when integrated over the energy gives us the total 
persistent current. Experimental detection of this persistent current would be via a measurement of the magnetic 
moment of the ring. 

Finally, we consider the case of a similar setting (as in Fig. 1(a)) but without any MBS. This system is a normal 
metal AB ring with an s-wave superconductor in its upper arm. The left lead is at potential Vi while no voltage is 
applied to the right. In Fig. 6 we plot the results for the non-local conductance G„; = (e^/7i)[l — i?^ + i?^ + Te ~ T^]- 
To conserve currents on either side of the superconductor one lets the potential in the superconductor float. Due to 
the absence of MBS the nonlocal conductance is symmetric with respect to magnetic fleld reversal (see top panel of 
Fig. 6). In the bottom panel of Fig. 6, we see Gni{Vi) 7^ Gni{~Vi). This is in contrast to the case where MBS are 
present. 



CONCLUSIONS 



We have introduced a novel mechanism to detect MBS occurring in TIs. The reahzation and control of MBS would 
be the first step towards a fault tolerant quantum computer. In Table I we summarize our results. We see that 
coupled MBS break magnetic flux symmetry. One also sees that period doubling occurs in presence of an electric 
potential and in the presence of individual MBS (Fig. 3). These can be easily used as means to detect MBS. Further 
there is a pronounced zcro-cncrgy dip/crest in presence of individual MBS which changes into a continuous function 
(without a maximum or minimum) in their absence. Finally, coupled MBS, can induce a persistent current in absence 
of a magnetic flux in a topologically insulating loop. Although these results are derived for a loop made of a TI, 
the results would remain valid in case the TI is replaced with a semiconductor provided a MBS can exist in such 
a system^U, since we used a generic Hamiltonian for the Majorana scatterer, see Eq. ([2|). An extension of the work 
would be to study the shot noise generated in our settings^. The symmetries of this could also be a matter of interest 
in the process of detection. 
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